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Abstract 

It is shown that when higher-shell admixtures are included for sys- 
tems with two valence particles or holes, there are effects which are 
quite different from those for one-valence-nucleon systems. For exam- 
ple, for nuclei with one valence particle or hole, there is no first-order 
correction for the magnetic dipole moment or the Gamow- Teller transi- 
tion amplitude. However for nuclei with two valence particles or holes, 
one can get substantial corrections. The effects of the tensor and spin- 
orbit interactions in core renormalization are emphasized. We find that 
in 6 Li, the spin-orbit interaction causes the quadrupole moment of the 
J 7r =l + ground state to be positive, but the tensor interaction causes 
it to be negative. The G-matrices derived from realistic interactions 
are employed in these calculations. 



1 Introduction 



There has been considerable study of core renormalization for a system of 
a closed LS shell plus or minus one nucleon. The most studied problem 
is probably the electric quadrupole (E2) effective charge. When doing cal- 
culations of E2 properties, i.e., quadrupole moments or £?(E2)'s in a single 
major shell, it is necessary to use effective charges for the valence protons 
and neutrons. The quadrupole moment for the ground state of 17 is about 
-3e fm 2 , suggesting that the effective charge for the valence d 5 / 2 -neutron is 
about -3/(-6)=0.5, where the quadrupole moment for a d5/2-proton is taken 
to be -6e fm 2 . The same analysis would lead to an effective charge of 1.2 for 
a valence proton, although more careful considerations may lead to some- 
what different values. The use of effective charges can be partly justified 
in first-order perturbation theory by the core-polarization process shown in 
Fig. 1(a). Larger effective charges can be obtained in R.P.A. calculations O. 
A typical RPA diagram is shown in Fig. 1(b). Although there have been cal- 
culations with certain interactions which show quantitative agreement with 
empirical values, there have been complete second-order perturbation calcu- 
lations by Siegel and Zamick O and Ellis and Siegel ]2[ which yielded very 
small polarization charges. Evidently the non-RPA diagrams cancel with the 
RPA diagrams. More recently, Skouras and Miither also obtained small 
polarization charges in higher-order calculations with G-matrices evaluated 
from Bonn potentials [§]]. 

At the extreme it is known that for a closed LS shell plus or minus 
one particle, the first-order perturbation-theory corrections for the magnetic 
dipole (Ml) operator and to the Gamow- Teller (GT) operator [\J~^{gil+g s s) 
and ot± respectively] vanish. One has to go to higher-order perturbation 
theory and invoke exchange currents in order to explain deviations from the 
lowest order. 

In this work, we will study higher-shell corrections to properties of sys- 
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terns involving one and two valence nucleons or holes. Although the first- 
order perturbation-theory correction vanishes for systems with one valence 
particle or hole, we may get some effects for systems with two valence parti- 
cles or holes. For example, the effective interaction between the two valence 
nucleons or holes can be renormalized (see, for example, Fig. 2) and hence 
the two-body wave function can change. 

Rather than perform perturbative calculations, we will perform matrix 
diagonalizations using the OXBASH shell-model (SM) code |^]. The model 
space will consist of one or two valence particles or holes, plus excitations 
of the core through "27ia;" . By this we mean that either two core- nucleons 
will be excited through one major shell or one core-nucleon will be excited 
through two major shells. 

There have been previous SM calculations with "2hu n corrections. For 
example, van Hees et al || and Wolters et al jjj have extended the clas- 
sic "Ohuj" calculations of Cohen and Kurath Q by including higher shells. 
At the time of this writing, another paper has appeared |J. Whereas the 
above works have used a phenomenological interaction and focused on get- 
ting a global fit to experiment — certainly a worthwhile quest — we are 
more concerned with seeing how well the G matrices, derived from modern 
realistic nucleon-nucleon (NN) potentials, are able to account for the mea- 
sured properties. Furthermore, we want to see how the various parts of the 
NN interaction, especially the tensor and spin-orbit components, affect the 
physical observables. 

The calculations are performed with Brueckner G-matrices |l0| calcu- 
lated from a new Nijmegen potential (Nijmll) JllJ . Comparisons with other 
potentials, including the Hamada- Johnson hard-core ]l2|, the original Reid- 
soft-core and a new Reid-soft-core (Reid93) |ll|] potential, will be made 
on a selected basis. It should be emphasized that our calculations will not 
assume an inert core. The single-particle energies are implicitly generated in 
OXBASH H from the two-body G-matrix elements as well as the one-body 
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kinetic energy. More explicitly, the matrix diagonalization will be performed 
for the SM Hamiltonian 



(A \ A ,, 

H SU = [J2 ** - T CM + I] Gij + \{H CM ~ ^hu), (!) 
Vi=l / i<i 



where t and Tcm are the one-body and center-of-mass (CM) kinetic energies 
respectively. The G matrix is calculated according to 

G{E S ) = vi 2 + v l2 - Q rv n , (2) 

E s - (/ii + h 2 + V12) 

with £J S the starting energy and hi = U + Ui = + |mu; 2 r? the one-body 
harmonic-oscillator Hamiltonian. We choose ?jo;=16MeV for A=f> and 6 
and hu=\AM(N for A=1A, 15, 17 and 18. The Pauli operator Q in Eq.® is 
defined to exclude the scattering of the two particles in the ladder diagrams 
into the occupied states as well as the states that will be included in the 
model space. For example, in the ladder diagrams, the scattering into an 
intermediate state with two nucleons in the ls-Od major shells or a state 
with one nucleon in the Op shell and the other in the Of-lp shell is forbidden 
(Q=0), since these states will be included in our (0 + 2)hu shell model 
diagonalization. The starting energy E s in the G-matrix (|2|) is taken to be 

E s = e a + e b + A, (3) 

where (e a + e&) is the unperturbed energy for the initial two-particle state 
in the ladder diagrams [e» = (2nj + U + \)huj with i = a, b] and A, roughly 
speaking, represents the interaction energy which in this work is treated as 
an adjustable parameter whose value is chosen to yield a reasonable nuclear 
binding energy. 

Note that the G-matrix elements used in the Ohoj calculation will be the 
same as those used in the (0 + 2)htu calculation so that the difference in the 
results is solely due to the core-renormalization effects. 
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We force, for the lower-lying states, the spurious CM motion to be in its 
lowest Ofuj configuration by adding to the SM Hamiltonian the last term in 
Eq.© with A > 1. 



2 Results for systems with one valence particle or 
hole 

We give results using the G-matrices from the Nijmegen NN potential j[l]] 
for A=5, 6, 14, 15, 17 and 18 in Table I. The data in moments is taken from 



the compilations of Ajzenberg-Selove [14] and Raghavan 



We emphasize again that the single-particle energies in our SM studies 
are not taken as parameters, rather, they are implicitly generated by the 
kinetic energy and the G-matrix. In a (0+2)fiuj space, the calculated "single- 
particle" splitting between the 1/27/ and 3/2j~ states in 5 He is 3.633 MeV. 
In 15 0, the 3/27/-1/27/ splitting is 5.479 MeV. The splitting between the 
3/2+ and 5/2+ states in 17 is 6.260 MeV. 

Focusing first on the one-valence-nucleon system A=17, we see that in- 
deed the values of B(GT) and B(M1) are nearly the same in the Ohuj space, 
in which the configuration is a d$/2 nucleon outside a closed 16 core, as in 
the larger (0 + 2)htv space. The value of B(GT) [our definition is such that 
the factor (1.251) 2 is not included; the operator is ^ °"(*)^±(*)] i n the small 
0hcu space is 1.400 and in the large (0 + 2)hu> space 1.381. For ^4=15, the 
story is the same. The value of -B(GT) changes very little (from 0.333 to 
0.326) in going from the Otku to the (0 + 2)Huj space. Underlying this be- 
havior for ^4=15 and A=17 is the previously mentioned theorem that there 
are no first-order corrections to the GT or Ml moment for a closed LS shell 
plus or minus one particle. The small deviations are due to the fact that 
the shell model goes beyond first order. 

We next consider the E2 properties for A=17. In the simple SM picture, 
the quadrupole moment of 17 is zero because the valence nucleon is a 



5 



neutron. In first-order perturbation theory, there is a contribution due to 
the excitation of a core proton through two major shells, e.g. from Os to 
ls-Od or from Op to Of-lp. In the SM calculation that we do here, that effect 
is implicitly taken into account, but also one has the excitation of two core 
nucleons through one major shell, i.e., (0p)~ 2 (ls0d) 2 . These configurations 
contribute only in second-order perturbation theory. 

The calculated polarization charge for the d5/2-neutron is the ratio of 
the quadrupole moment of 17 to that of a G?5/ 2 -proton which can be read 
off from Table I: 

-0.720 

0.123. (4) 



-5.924 

Similarly, for the valence proton, 



-6.183 - (-5.924) _ ^ 



-5.924 

These are much smaller than the phenomenological values of 0.5 and 0.2. 
These results agree qualitatively with the higher-order perturbation theory 
results in Ref.|| where small polarization charges are also obtained. 

To check if the second-order contributions are canceling the first-order 
ones, we consider the first-order perturbation theory. If we use an energy 
denominator of (—2hu) for Fig. 1(a), the polarization charge for the neutron 
is e v = 0.385 and that for the proton is = 0.131. However, when the 
calculated single-particle energy splittings are used, the values are e u = 0.263 
and e n = 0.087. These are still larger than the SM values of 0.123 and 0.044. 

These discrepancies are partly due to the fact that in the above first- 
order perturbation-theory approach, the wave function is not normalized. 
The ground-state wave function of 6 Li in the (0 + 2)hu SM calculation is 

|g.s.) = 68.446% \lp-0h) + 27.727%|2p- Ih) + 3.827%|3p-2/i). (6) 

A properly normalized wave function in the perturbation theory would there- 
fore lead to an approximately 32%'s reduction (from 0.263 to 0.18) in the 
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first-order contribution to the effective charge. Further reductions are caused 
by the higher-order effects that are included in the (0 + 2)Huj SM calculation. 

We thus see that one must exhibit care when employing the perturbation- 
theory approach to calculate effective operators. The first-order perturba- 
tion theory first appears to be accountable for the large empirical polariza- 
tion charge of nearly 0.5 for the d5/2-neutron in 17 0, a more careful analysis 
involving the use of the self-consistent single-particle energies and the nor- 
malized wave function results in a factor of two's reduction in e v , leading 
to a much smaller result (0.385 to 0.18). The (0 + 2)hu SM value of 0.123 
is more consistent with the refined first-order perturbation-theory result of 
0.18, both are smaller than the empirical value. The large quadrupole mo- 
ment observed for 17 is due to the fact the "spherical core" is actually 
deformed and carries a non-zero angular momentum. To give sufficiently 
large polarization charges, bigger than (0 + 2)huj calculations are needed 
which will simulate highly deformed admixtures into the ground state. 

3 Higher-shell effects for systems with two valence 
nucleons (or holes) 

For A=1A, we have a striking example of where there is a large effect 

due to higher-shell admixtures, while there is almost none for a one-valence- 
nucleon system. The value of B(GT) for the decay 14 C( J = + , T = 1) — > 14 
N( J = 1+, T = 0) changes from 2.627 to 1.305. The value of B(M1) changes 
by almost the same ratio; indeed from isospin consideration the spin part 
of £>(M1) and B(GT) must change by exactly the same ratio. The large 
change has been previously noted by us [|16| . 

What must be happening is that the wave function of the two valence 
nucleons is changing due to a renormalization of the effective interactions 
from higher shells. A typical second-order diagram is shown in Fig. 2. We 
recall that this weak beta decay is famous because, although it is a prior an 
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allowed GT transition, the life-time for the transition is very long, indicating 
a near vanishing GT matrix element. 

Furthermore, it was noted by Inglis [17| that in a two-hole calculation 
one needs a tensor interaction to get a vanishing GT matrix element. In 
our calculation we obtain a B(GT) value of 1.305. We believe that the 
bare tensor interaction is too strong so if we imagine turning on the tensor 
interaction from zero to its full (bare) value, we get an overshoot. The matrix 
element first gets smaller in magnitude, then goes through zero, changes sign 
and gets large in magnitude again. The effect of the higher-shell admixtures 
is to effectively weaken the tensor interaction in the valence space, thus 
causing the GT matrix element to decrease in magnitude. 

For A=6, there are some encouraging results from higher-shell effects. 
The experimental magnetic dipole moment and electric quadrupole moment 
of the 1 + ground state of 6 Li are respectively 0.822/^ and -0.082efm 2 



[18|. The sign of the 6 Li quadrupole moment is opposite to that of the 
deuteron (Q = +0.28efm 2 ). In the LS limit, the ground state of 6 Li has 
a p 2 configuration with L=0, 5=1, J=l. In that limit, the Ml moment is 
M = + fJ-p = (—1.913 + 2.793)/z/v = 0.880/iAr and the quadrupole moment 
is zero. 

In the small space, the results for fi and Q are respectively 0.866/Utv and 
-0.237efm 2 . In the large space, the results are 0.849^at and -0.163efm 2 in 
both cases closer to experiment. 

In Table II we give the results of the jj and LS limits for the Ml and 
E2 moments of the 1 + ground state of 6 Li. For the Ml moment, these two 
limits yield respectively = (^p + ^ n+1 ) = 0.627/uat and fi = [i p + fi n = 
0.880/ijv. The experimental value 0.822/xtv lies between these two extremes. 
The quadrupole moment Q is zero in the LS limit (L=0, S=l, J=l) and 
is positive in the jj limit. The expression for Q in the jj limit can be 
related to that of a p 3 / 2 proton via the Wigner-Eckart theorem in terms of 
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Clebsh-Gordan and unitary Racah coefficients: 

= -0AQ(p 3/2tn ). (7) 

The quadrupole moment of a j = I + 1/2 proton is —lh/(mcu). Hence the 
quadrupole moment for a p 3 / 2 -proton is negative: Q(p3/2,n) = —h/imuj). 
The (p^/2) J=1 ' T= ° quadrupole moment is therefore positive (1.04 efm 2 when 
fiuj = 16 MeV). 

To get further insight into the A=6 behavior, we use a schematic inter- 
action 

^sche = V C + xV so + yV t . (8) 

With x=y=l, this interaction gives relative matrix elements similar to those 
of Bonn A Q for the region of 16 UlSfl . Only partial waves with a relative 
angular momentum /=0 and 1=1 for the diagonal and /=0 to 1=2 for the 
off-diagonal are included |l6|. It is constructed so that it is easy to control 
the strengths of the spin-orbit and tensor interactions by adjusting x and 
y, e.g., for x=0, we turn off the spin-orbit interaction and for y=0, we turn 
off the tensor interaction. The results for the schematic interaction are also 
presented in Table II. The matrix elements of the interaction have been 
multiplied by 1.15 since they were fitted to 16 but are being applied to 
A=6. 

We verify that for x=y=0, one gets the LS wave function with L=0, 
S=l for the ground state of 6 Li. The magnetic moment in both the OHuj 
and (0 + 2)hu> calculations is \i n + /i p =0.880^Ar and the quadrupole moment 
is zero. These results transcend the truncated shell model. 

Note that the spin-orbit interaction and the tensor interaction have op- 
posite effects on the quadrupole moment. In the small space, when the spin- 
orbit interaction is turned on at full strength and the tensor interaction is 
turned off (i.e., x=l, y=0), the quadrupole moment Q is 0.121efm 2 . This 
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is consistent with the fact that making the spin-orbit interaction stronger 
takes one towards the jj limit of Table II. 

In the opposite extreme with the tensor interaction on and the spin-orbit 
interaction off (i.e., x=0, y=l), the value of Q is -0.608efm 2 . For x=y=l, 
the value of Q is -0.281efm 2 . In a previous work |2C], we argued that the 
effective spin-orbit interaction in the p shell should be even stronger than 
that given by the relativistic Bonn A with m*/m=l Q. A reasonable set 
of parameters is x=1.5, y=l with which the small-space value of Q is - 
0.114efm 2 , in reasonable agreement with experiment. The main point here, 
however, is that the quadrupole moment of 6 Li is extremely sensitive to the 
relative strengths of the spin-orbit and tensor interactions in the nucleus. 

As seen in Table I, the 2%tu admixtures enhance the quadrupole proper- 
ties of A=5, as expected. But but for ^4=6, the magnitude of the quadrupole 
moment decreases from -0.255 to -0.175efm 2 when 2hco configurations are 
included. We believe that for A=6 there are two opposing tendencies. The 
building up of the effective charge would cause the magnitude of Q to in- 
crease. But the "self weakening" of the tensor interaction due to higher-order 
admixtures will cause Q to go from negative towards positive values. 

A cluster approach for 6 Li, i.e., an alpha-deuteron cluster model calcu- 



lation has been carried out by Eskondarian et al [21]. Without going into 
too much detail we focus on their results. Whereas we get the magnetic 
moment of 6 Li larger than experiment, they get it smaller. Whereas we get 
the quadrupole moment of 6 Li to be negative, they get it positive. It is 
not clear if the differences are due to the differences between the SM and 
the cluster model, or in the differences between the interactions that are 
used. The SM will also yield a positive quadrupole moment if the spin-orbit 
interaction is stronger and/or the tensor interaction is weaker than what 
we have. Another possibility to be explored is that the cluster model has 
correlations which require even larger SM calculations than the ll 2hoj" ones 
performed here. 



10 



4 Closing remarks 



We have shown not only that higher-shell effects are important for nuclear 
properties but that they enter in a nonuniform way as we add nucleons to 
the nucleus. For example, higher-shell effects are much larger for the spin 
properties of a two-valence-particle (or hole) system than they are for a one- 
valence-particle (or hole) system. As we go from a small space (Ohu) to a 
large space [(0+2)7lu/|, the value of B(GT) for ^4=15 experiences little change 
(from 0.333 to 0.326), but for A=1A the change is much bigger (from 2.627 
to 1.305). This can be understood by the fact that the effective interaction 
between the two valence holes gets renormalized. In this renormalization, 
the effective tensor interaction (acting in the valence space) becomes weaker. 
The A=6 and A=14 cases offer good tests of the effective interaction. For 



6 Li, which has often been proposed as an ideal isoscalar target (see Ref. [18] 
for more details), the results for the J=l + static quadrupole moment and 
magnetic moment are very sensitive to the details of the spin-orbit and 
tensor interactions that are used. There is a surprisingly large scatter of 
theoretical results — both positive and negative values for the quadrupole 
moment of the J=l + ground state of 6 Li have been obtained. We note that 
the tensor interaction gives a negative value and the spin-orbit interaction 
a positive value for this quantity. Our results with the effective interactions 
calculated from the Nijmegen potential Nijmll as applied to "two-particle" 
systems are encouraging in the sense that higher-shell effects consistently 
go in the right direction. There is still a problem, however, in that the 
one-body quadrupole polarization charges come out too small. 

If we want to use the nucleus as a laboratory for fundamental studies, 
e.g., of fundamental symmetries, applications to astrophysics, or, as we have 
emphasized here, the effective interaction between two nucleons in a nucleus, 
we must include higher-shell effects. Indeed, our work here represents only 
a step in the right direction. In the future, excitations beyond '2/iu;' will 
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have to be included. Still, it is encouraging that we are able to get a better 
systematic feeling for how higher-shell admixtures affect nuclear properties 
in a nuclear medium. 
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Table I. The static and transitional properties in selected nuclei with one or two 
valence nucleons (or holes) in the Ohuj and (0 + 2)?ia; shell-model calculations using 
G-matrices calculated from a new Nijmegen local NN potential pj . B(M1) is in 
units of fj/ff, fx in units of /x^r and Q in units of efm 2 . Experimental data are taken 
from Refs.Q ||. 



Nucleus 


Quantity 


Ohuj (0 + 2)?kJ 


Expt 


A=5 


B(GT)(3/2- -> 3/2") 

^( 5 Hc) 

M( 5 Li) 

Q( 5 He) 

Q( 5 Li) 


1.666 1.632 
-1.913 -1.864 
3.793 3.741 
-2.592 -2.722 
0.000 -0.303 




A=6 


£(GT)(0+1 -» 1+0) 
£(M1)(0+1 -► 1+0) 
M (l+0) 
Q(l+0) 


5.626 5.397 
16.29 15.57 
0.866 0.849 
-0.255 -0.175 


4.847 

0.822 
-0.082 


A=U 


B(GT)(0+1 -> 1+0) 
B(M1)(0+1 -> 1+0) 
M (l+0) 
Q(l+0) 


2.627 1.305 
6.916 3.842 
0.634 0.499 
1.761 2.164 


0.000 

0.404 
1.56 


A=15 


B(GT)(l/2~ -» 1/2") 
^( 15 N) 

m( 15 o) 


0.333 0.326 
-0.264 -0.276 
0.638 0.654 


0.270 
-0.283 
0.719 


A=n 


B(GT)(5/2+ -» 5/2+) 

m( 17 o) 
m( 17 f) 
Q( 17 o) 

Q( 17 F) 


1.400 1.381 
-1.913 -1.869 
4.793 4.748 
0.000 -0.720 
-5.924 -6.183 


1.089 
-1.894 
4.722 
-2.578 
-5.8 


A=18 


B(GT)(0+1 -> 1+0) 
B(M1)(0+1 -► 1+0) 
M (l+0) 
Q(l+0) 


5.312 5.436 
16.75 16.19 
0.858 
-0.514 


3.301 
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Table II. The magnetic dipole moment /i (in units of hn) and quadrupole moment 
Q (in units of efm 2 ) for the ground state of 6 Li in the Oftw and (0 + 2)huj shell- 
model calculations using (1) the schematic interaction with varying spin-orbit (x) 
and tensor (y) strengths, and (2) G-matrices calculated from different realistic NN 
potentials: Hamada-Johnson jL2|, Reid-soft-core new Reid-soft-core (Reid93) 
[ [il"|| , and new Nijmegen fli"|| . 



Interaction 


A* 


Q 


X 


V 


ftT>t,i 






(ft 4- 9\f>i,t 


0.0 


0.0 


u.oou 


u.oou 


n ono 

u.uuu 


u.uuu 


0.0 


0.5 


0.867 


0.859 


-0.345 


-0.415 


0.0 


1.0 


0.840 


0.816 


-0.608 


-0.711 


0.5 


0.0 


0.877 


0.878 


+0.038 


+0.033 


0.5 


0.5 


0.872 


0.866 


-0.221 


-0.276 


0.5 


1.0 


0.856 


0.835 


-0.446 


-0.533 


1.0 


0.0 


0.869 


0.871 


+0.121 


+0.106 


1.0 


0.5 


0.870 


0.866 


-0.085 


-0.135 


1.0 


1.0 


0.863 


0.845 


-0.281 


-0.358 


1.5 


0.0 


0.854 


0.859 


+0.221 


+0.195 


1.5 


0.5 


0.861 


0.860 


+0.055 


+0.006 


1.5 


1.0 


0.860 


0.848 


-0.114 


-0.185 




LS limit 




0.880 




0.000 




33 limit 




0.627 




1.037 


Hamada-Johnson 


0.861 


0.840 


-0.364 


-0.285 
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Figure caption 



Fig.l Core-polarization diagrams for E2 effective charge, (a) First-order 
diagram; (b) Typical RPA graph. 

Fig. 2 Core renormalization of the effective interaction. This changes the 
two-particle wavefunction and hence static properties like Ml and E2 mo- 
ments. 



17 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/93 1 1016vl 



